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^ ■ ABSTRACT 

^ ' Spectrophotometric observations of the gravitationally microlensed quasar 

Q2237+0305 during a High Magnification Event (HME) is potentiaUy a very power- 
ful tool for probing the structure of the quasars accretion disc on scales of ~ 10~^ 
arc seconds. In cases where the HME is produced by a single caustic (SHME), mi- 
crolensing induced changes in the spectrum during the event may be used to directly 
' infer the source intensity profile (eg. Grieger et al. 1988; Yonehara et al. 1988; Agol 

^ \ & Krolik 1999a). Several groups are actively monitoring Q2237-f0305 with this goal 

^— ^ ■ in mind. How often we can expect to observe a HME is dependent on the lens system 

^ ^ ' parameters of galactic transverse velocity, mean microlens mass and the size of the 

magnified continuum source. We have previously used published microlensed light- 
' curves to obtain expressions for the likely-hood of the values for these parameters 

Qs, \ (Wyithe, Webster & Turner 1999b,c; Wyithe, Webster, Turner & Mortlock 1999). 

Qv^ . Here we use this information to investigate the expected rate of SHMEs. We find the 

average rate of SHMEs as well as the number that we can expect to observe over 
periods of a decade and of a single observing season. We find that the average SHME 
rate summed over all images in Q2237-I-0305 is 1.5±0.6 - 6.2±1.3 events per decade. 
O ' During the period following a caustic crossing we find that the event rate in the corre- 

|H \ spending image is enhanced by ^ 50 — 100%, and therefore that the overall event rate 

c/3 ■ may be higher during these periods. From the distribution of events expected during 

a 6 month period we find that there is 1 chance in 4 — 10 of observing a SHME per 
observing season. The systematic dependence in these values arises from the different 
assumptions for smooth matter content, orientation of the galactic transverse velocity 
. ^ and the size of photometric error in the monitoring data. The results support contin- 

^ . ued monitoring of Q2237+0305 with the aim of obtaining detailed spectroscopic and 

photometric observations of a SHME. 

Key words: gravitational lensing - microlensing - numerical methods. 



1 INTRODUCTION 

The object Q2237+0305 (Huchra et al. 1985) comprises a 
source quasar at a redshift oi z = 1.695 that is gravitation- 
ally lensed by a foreground galaxy with z — 0.0394 pro- 
ducing 4 resolvable images with separations of 1". Each 
of the 4 images are observed through the galactic bulge, 
which has a microlensing optical depth in stars that is of 
order unity (eg. Kent & Falco 1988; Schneider et al. 1988; 
Schmidt, Webster & Lewis 1998). In addition, the proximity 
of the lensing galaxy means that the effective transverse ve- 
locity may be high. The combination of these considerations 
make Q2237-I-0305 the ideal object from which to study mi- 
crolensing. Indeed, Q2237-(-0305 is the only object in which 



cosmological microlensing has been confirmed (Irwin et.al 
1989; Corrigan et.al 1991). 

During a microlensed high magnification event (HME), 
regions of the source that are much smaller than the mi- 
crolens Einstein radius are differentially amplified. Simula- 
tions by Wambsganss & Paczynski (1991) showed that this 
produces an observed colour change that should be observ- 
able. Irwin et al. (1989) found evidence for such a colour 
change during the 1988 event in image A of Q2237-I-0305. 
Later Lewis et al. (1998) determined the ratios of emission 
line equivalent widths relative to one image. They show that 
the ratios vary between images at a single epoch, and that 
the ratio for a single image (image A) varies between two 
different epochs. This spectral change is interpreted as be- 



© 1999 RAS 



2 J. S. B. Wyithe et al. 



ing due to the different spatial extents of the continuum 
and emission hue regions being differentially amplified due 
to microlensing. 

While reverberation mapping provides a technique for 
determining the geometric structure of a QSO, typical 
sampling time-scales are several days which correspond to 
2.5 X 10^^ cm/day, comparable to the expected size of the ac- 
cretion disc. Observation of a SHME provides a probe of the 
quasars central engine on scales at least an order of magni- 
tude smaller than this. Two approaches have been taken to 
study microlensed quasar spectra. Firstly, several authors 
have calculated the effect of microlensing on the spectra 
and luminosity of model continuum regions (eg. Jaroszynski 
et al. (1992); Ranch & Blandford (1992); Yonehara et al. 
(1998)). A different approach is to invert the light curve ob- 
served to directly obtain a surface brightness profile of the 
source (eg. Grieger et al. 1998, 1991; Agol & Krolik 1999). 
Which ever approach is taken, good time resolution of ob- 
servations (~ 1 per day), and broad wavelength coverage are 
required for the duration of the SHME so that the contin- 
uum region can be probed with the finest possible resolution 
over a large range of spatial scales. It is therefore essential to 
frequently monitor the quasar in order to detect the onset of 
the HME as early as possible. Daily ground based monitor- 
ing of Q2237-I-0305 is underway at the Apache point obser- 
vatory. In addition monitoring is implemented by the OGLE 
collaboration. The aim of the former program is to observe 
the onset of an HME so as to obtain spectra at different 
times during the event. 

The rate of HMEs has been investigated by several 
authors (eg. Wambsganss, Paczynski & Schneider (1990); 
Wambsganss, Paczynski & Katz (1989); Witt, Kayser & 
Refsdal 1993). However the expected rate of events has been 
poorly constrained due to the absence of information on sev- 
eral microlensing parameters. In particular the event rate 
is proportional to the effective galactic transverse velocity, 
and inversely proportional to the square root of the mean 
microlens mass. The size of the source will also effect the 
proportion of HMEs that are SHMEs. In this work we cal- 
culate HME and SHME rates at many assumed values for 
the above parameters, and convolve these rates with proba- 
bilities for the parameter values (Wyithe, Webster & Turner 
1999b,c; Wyithe, Webster, Turner & Mortlock 1999). 

This paper is presented in four parts Sections ^ and ^ 
describe the microlensing model for Q2237-f 0305 and the 
calculation of the expected event rates respectively. Section 
k describes the results obtained. 



2 THE MICROLENSING MODEL 

Throughout the paper, standard notation for gravitational 
lensing is used. The Einstein radius of a IMq star in 
the source plane is denoted by 770- The normalised shear 
is denoted by 7, and the convergence or optical depth 
by K. Kc and k, are the optical depth in smoothly dis- 
tributed matter and stars respectively. To construct a mi- 
crolensed light-curve we use the inversion technique of Lewis 
et al. (1993) and Witt (1993). For the microlensing mod- 
els of Q2237-I-0305 presented in the current work we as- 
sume the macro-parameters for the lensing galaxy calcu- 
lated by Schmidt, Webster & Lewis (1998). These values 



Table 1. Values of the total optical depth and the magnitude of 
the shear at the position of each of the 4 images of Q2237-f 0305. 
The quoted values are those of Schmidt, Webster & Lewis (1998). 
K, and Kc are the optical depths in stars and in smoothly dis- 
tributed matter respectively. 



Image 




l7l 


A 


0.36 


0.40 


B 


0.36 


0.40 


C 


0.69 


0.71 


D 


0.59 


0.61 



are shown in table y. Where required a cosmology having 
SI = 1 with Ho = 75 km sec~^ is assumed. For our analysis 
we assume that microlensing is produced through the com- 
bination of a galactic transverse velocity with an isotropic 
velocity dispersion. Moreover, we assume that the magni- 
tude of the dispersion is the same at the position of each of 
the four images. For the line-of-sight velocity dispersion of 
the stars in the galactic bulge we take the theoretical value 
of CT* ~ 165 fern sec~^ (Schmidt, Webster & Lewis 1998). 

We describe the microlensing rate in terms of the effec- 
tive transverse velocity which is defined as the transverse ve- 
locity that produces a microlensing rate from a static model 
equal to that of the observed light-curve (Wyithe, Webster 
& Turner 1999b). The effective transverse velocity therefore 
describes the microlensing rate due to the combination of 
the effects of a galactic transverse velocity and proper mo- 
tion of microlenses. When calculating HME rates we assume 
that the effective transverse velocity accurately describes not 
only the time averaged HME rate, but also the temporal 
clustering of HMEs and the distribution of orientations be- 
tween the caustic and source trajectory. 

We have previously obtained the following normalised 
probability distributions. 

i) Pa{Ss\{m),Veff), the probability that the continuum 
source diameter is between Ss and 5*5 + dSs given a mean 
microlens mass (m) and an effective galactic transverse ve- 
locity v^ff (Wyithe, Webster, Turner & Mortlock 2000). 

ii) Pv{veff\{ni)) the probability that the effective galactic 
transverse velocity is between Vef / and Vsf / -I- dvef / given 
a mean microlens mass (m) (Wyithe, Webster & Turner 
1999b). 

iii) Pmiim)), the probability that the mean microlens mass 
is between (m) and (m) + d(m) (Wyithe, Webster & Turner 
2000). 

Pv{vef f\{m)) and pm{{m}) were computed using flat 

(piVtran) <X dVtran), aud logarlthmlc (piVtran) OC ^^^) 

assumptions for the Bayesian prior for galactic transverse 
velocity (Vtran). Pv{vef f\{m)) was found to be insensitive to 
the prior assumed, however pm((m)) showed some depen- 
dence. For a point source, the rate of events calculated from 
Pviveff\{m)) is independent of Pmi{m)), however a slight 
dependence is introduced due to the finite size of the source 
in relation to the caustic network. In the remainder of this 
paper we use pmiim)) as calculated using the assumption of 
a logarithmic prior. We note that the assumption of the flat 
prior raises the computed event rates by a few percent. 

The functions Ps{Ss\{m),Veff), pv{veff\{ni)), Pm{{m)) 
and the event rates presented in this paper were computed 
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for the following assumptions of smooth matter density, pho- 
tometric error, and direction of the galactic transverse ve- 
locity. Two models are considered for the distribution of mi- 
crolenses, one with no continuously distributed matter, and 
one where smooth matter contributes 50% of the surface 
mass density. Two orientations were chosen for the trans- 
verse velocity with respect to the galaxy, with the source 
trajectory being parallel to the A— B and C— D axes. The 
two orientations bracket the range of possibilities, and be- 
cause the images are positioned approximately orthogonally 
with respect to the galactic centre correspond to shear values 
of 7A, 7b < 0, 7c, 7i3 > and 7a, 7b < 0, 7c, 7d > respec- 
tively. The simulations used two different estimates of the 
error in the photometric magnitudes. In the first case a small 
error was assumed (SE). For images A and B, crsB=0.01 
mag, and for images C and D asE=0.02 mag. In the second 
case, a larger error was assumed (LE). For images A and B, 
aLE=0-02 mag and for images C and D (tlb=0.04 mag. The 
observational error in Irwin et al. (1989) was 0.02 mag. 

Both the microlensing rate due to a transverse veloc- 
ity (eg. Witt, Kaiser & Refsdal 1993), as well as the corre- 
sponding rate due to proper motions (Wyithe, Webster & 
Turner 1999a) are not functions of the details of the mi- 
crolens mass distribution, but rather are only dependent on 
the mean microlens mass. We therefore limit our attention 
to models in which all the point masses have the same mass 
since the results obtained will be applicable to other models 
with different forms for the mass function. The microlensing 
models presented in this work have been discussed in detail 
in Wyithe, Webster & Turner (2000). 

We obtain caustic positions and orientations as follows 
using the contouring algorithm (Witt 1993; Lewis et al. 
1993), which we implement according to the description of 
Lewis et al. (1993). The image curves are followed in discrete 
steps in one direction. During the contouring algorithm the 
situation of the image line crossing a critical curve results 
in the image magnification changing sign due to a parity 
flip between the critical images on either side of the critical 
curve. The location of the critical curve can then be found 
through a search for the image position that produces a lens 
mapping Jacobian determinant of zero. This value for the 
critical image position is used to calculate the corresponding 
caustic position. The gradient of the caustic tangent is found 
by first computing the location of a point a short distance 
along the critical curve in both directions. The correspond- 
ing caustic points are then used to compute the gradient of 
the caustic tangent from a three-point derivative. 



3 EVENT RATE 

We find the rate of caustic crossings as a function of mean 
microlens mass, source size and transverse velocity. We de- 
fine two classes of HMEs for this purpose: 1) those that 
involve a single caustic and so can be deconvolved to yield 
information on the source intensity profile etc. of the source 
(SHME), and 2) those where the source is directly effected 
by two or more caustics at once (MHME). The following 
algorithm is used to automatically classify a model HME 
into one of these classes. We define a caustic to be involved 
in an event if the closest edge of the source is within one 
source diameter of the caustic on the -I- side (source has as- 



sociated critical images), or in contact with the caustic on 
the - side (source has no associated critical images). It is 
computationally prohibitive to compute full light curves for 
extended sources of many sizes over a statistically represen- 
titive number of starfields. Therefore, we make the follow- 
ing approximation for the small sources (<C Irio) that we 
are considering here. We assume that the caustic is locally 
straight over a length that is larger than the source. If the 
caustic makes an angle with the source track of 9, then a 
circular source comes into contact with it at a distance of 
L = „ ^ a along the source line from the source centre. This 

2 sin ^ 

length can be calculated for the nearest caustic along the 
source line in both directions {Lujt, Lright)- The source is 
only effected by one caustic at a time if the distance between 
these caustics along the source line is less than Lieft + Lright 
for H — , 25* + Lie ft + Lright for — h and S -\- Lie ft + Lright 

for ++ or caustic pairs. If this condition is fulfilled for 

a source on both sides of a caustic then the source is said to 
undergo a SHME when crossing that caustic. 

The average rate of SHMEs for a given source size, mi- 
crolens mass and effective galactic transverse velocity is 



R{S, (m» 



(1) 



where Ro {S, (m)) is the rate for Veff {{m)) = 1. The proba- 
bility of R {S, (m)) lying between R and _R + A_R is therefore 



PR 



r{R\S, (m)) 



Pv Veff 



Ro (S, (m)) 



|(m) 



Ro{S, (m)) 



(2) 



From this we can evaluate the probability for the SHME 
rate given our previous calculations for the probability of 
microlens mass and source size: 



PRiR) 



dm / dS Ps{S\{m),Veff)p,n{{m))pii(^R\S,{m)'j (3) 



where 



'Veff 



Ro {S, (m)) 



(4) 



The caustic network is clustered on temporal scales of 
decades. Therefore, the average event rate may not be ob- 
servable for as long as a century (Wambsganss, Paczynski 
& Schneider 1990). For this reason the distribution of the 
number of SHMEs expected during a finite period is of more 
utility. We calculate the histogram for the number of SHMEs 
per period for a collection of assumed effective galactic trans- 
verse velocities, source sizes and mean microlens masses: 



hN iN\S, (rn),Veff) . 



(5) 



The probability for finding a rmmber of SHMEs in a period 
can then be obtained: 

hN (N) = 



dm J dS J dveff {ps{S\{m),Veff)pm{{m)) 

X Vv{Veff\{m))hN{N\S,{m),Veff)) (6) 

PflT. {Ri) and /ijv; [Ni) are evaluated for each of the four 
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Table 2. The means and variances of the probabihty distribution functions for HME and SHME rates (events/decade). 







1A,B > 


7C,D < 






7A,s < 0,7c,D > 




Error (Am = ±0.02) 


0%Kc 


50%Kc 


0%Kc 


50%Kc 




HME rate 


SHME rate 


HME rate 


SHME rate 


HME rate 


SHME rate 


HME rate 


SHME rate 




7.3±1.5 


6.2±1.3 


5.6±1.3 


4.8±1.3 


4.3±1.2 


3.5±1.0 


2.8±0.8 


2.3±0.7 




4.8±1.5 


4.1±.4 


3.7±1.3 


3.2±1.1 


2.8±1.0 


2.3±0.8 


1.8±0.7 


1.5±0.6 



Table 3. The means of the histograms h]v(Af) for SHME rates (events/decade) in the cases of events during randomly chosen 10 year 
periods (Uncorr. rate) and 10 year periods imediately following a caustic crossing (Corr. rate). 







7A,S > 0, 7C,D < 






7A,B < 0, 7C,D > 




Error (Am = ±0.02) 


0%Kc 




50%K 


c 


0%Kc 




50%K 


c 




Uncorr rate 


Corr rate 


Uncorr rate 


Corr rate 


Uncorr rate 


Corr rate 


Uncorr rate 


Corr rate 


Image A/B 


2o"A,s, lo-c,D 


1.5 


2.3 


1.2 


1.9 


1.0 


1.7 


0.6 


1.2 


IfA.s, \^C,D 


1.0 


1.9 


0.9 


1.6 


0.7 


1.5 


0.5 


1.1 


Image C 


2o"A,B, lo"c,i3 


1.4 


2.2 


0.9 


1.5 


0.9 


1.5 


0.5 


1.1 




1.0 


1.9 


0.7 


1.4 


0.6 


1.3 


0.4 


1.1 


Image D 


2o-A,B, lo"C,D 


2.3 


3.3 


2.0 


2.7 


1.3 


2.0 


1.0 


1.6 


lo'A.B: 50"C,D 


1.7 


2.6 


1.4 


2.2 


0.9 


1.7 


0.7 


1.5 



images (i = 1 — 4). The total rate summed over all images 
is then evaluated as 

PRtot {Rtot) = J nUdR^pn^ (i?.) s {Rtot - {^UlR^)) .(7) 

Similarly, the probability for finding a number of SHMEs in 
a finite period summed over all images is 

^ivtot [Ntot) = J2 ntiftiv. (Nr) S {Ntot, StiiV,) . (8) 

Monitoring data (Lens Monitoring Project, Apache 
Point Observatory; OGLE collaboration ( see 



http: / / www.astro.princeton.edu / ogle /ogie2 /huchra.html ) ) 



shows evidence for recent high magnification events in two 
images. Clustering of caustics will tend to produce a higher 
than average rate of caustic crossings (in a single image) 
in the period following a HME. In addition to calculating 
/ijv (N) for randomly chosen 10 year periods, we therefore 
also calculate /ijv (A'^) for a sample 10 year periods that im- 
mediately follow a caustic crossing. 



4 RESULTS 

Table ^ contains the expectation values and variances for 
the total HME and SHME rate. The variances in individual 
distributions are comparable to the systematic uncertainty 
in the expectation value introduced by the assumptions of 
trajectory direction, photometric error, and smooth mat- 
ter content. The rates for HMEs and SHMEs are similar 
in all cases, with SHMEs being slightly less common. The 



small difference between these values is a result of the low 
probability for a source size that is > 0.1-y/ {m)r)o- The val- 
ues range between 1.5±0.6 and 6.2±1.3 events per decade. 
These rates are consistent with those calculated by Witt, 
Kayser & Refsdal (1993) who assumed a source plane trans- 
verse velocity of 6000 km sec~^ and a mean microlens mass 
of ~ O.35M0. The rates are lower where the trajectory is 
aligned with the A-B axis. This effect is larger than the low- 
ering of the event rate that is produced by the introduction 
of smooth matter. In addition, the predicted event rate is de- 
pendent on the photometric error assumed. This dependence 
comes from the systematic uncertainty that the assumption 
introduces into the determination of the effective transverse 
velocity. 

The microlensing parameter distributions were calcu- 
lated from an ensemble of 5000 simulations of the observed 
data for each image, and for each assumption of trajectory 
direction, photometric error, and smooth matter content. 
The 5000 simulations were spread over lOOOr^o of simulated 
light curves. If the error in the average rate is assumed to 
be due to Poison noise then it is less than ~ 1%. The er- 
rors in the modelling statistics are therefore not important 
in determining the probability distributions for the average 
HME rate. 

Figures ^ ^, p| and ^ display the functions pji{R), 
PRt^^iRtot), hN(N) and hNtot i^tot) obtained for the mi- 
crolensing models considered. Both the rates for the indi- 
vidual images as well as the combined rates are shown. 
fiN (N) and hNtot (Ntot) have been calculated for peri- 
ods of 1 decade, corresponding approximately to the total 
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Figure 1. The functions representing the probabiUty for event rates in Q2237+0305. Left: the time averaged rate, with expectation 
values, Centre: events expected per decade and Right: events expected per 6 month observing season for images A/B, C and D. The 
hght and dark lines correspond to 7 < and 7 > respectively, the thick and thin lines represent the resulting functions when the 
photometric error was assumed to be e=0-0^ mag in images A/B and (ts£;=0.02 in images C/D, and (t^s=0.02 mag in images A/B 
and (ri^B=0.04 in images C/D. The lens model contained no smoothly distributed matter. 




Figure 2. The functions representing the probability for combined image event rates in Q2237+0305. Left: the time averaged rate, with 
expectation values, Centre: events expected per decade and Right: events expected per 6 month observing season. The light and dark 
lines refer to assumed directions for transverse motion that are directed along the A-B axis and C-D axes respectively. The thick and 
thin lines represent the resulting functions when the photometric error was assumed to be (tsb=0.01 mag in images A/B and (Tg£;=0.02 
in images C/D, and crLE=0.02 mag in images A/B and crLE=0.04 in images C/D. The lens model contained no smoothly distributed 
matter. 



Q2237+0305 monitoring period, as well as a 6 month period 
which corresponds to a single observing season. 

We find that the event rates should be highest in im- 
age D, corresponding to the higher expected magnification. 
This is in contrast to the observed behaviour of image D 
which has been relatively faint for the period of monitoring, 
suggesting that the source may currently reside in a sparse 
region of the corresponding caustic network. It is important 
to note however that the non-observation of a HME in any 
image over the time that Q2237-I-0305 has been monitored 
should not be surprising, since the modes of all individual 
image histograms are zero. This offers an alternative to the 



explanation given by Witt & Mao (1994) that the lack of 
variability in images C and D (which persisted following the 
publication of that paper) is the result of an alignment be- 
tween the transverse velocity and the caustic clustering. 

The probability distributions of microlensing model pa- 
rameters that have been used to calculate the event rates 
have been obtained from the observed light curves. How- 
ever these statistics are derived primarily from the rate of 
low level microlensing variation which is related to but does 
not necessarily correspond directly with the rate of caustic 
crossings. It is therefore interesting to compare the number 
of HMEs in the observed light curves to the calculated distri- 
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Figure 3. The functions representing the probabiUty for event rates in Q2237+0305. Left: the time averaged rate, with expectation 
values, Centre: events expected per decade and Right: events expected per 6 month observing season for images A/B, C and D. The 
hght and dark lines correspond to 7 < and 7 > respectively, the thick and thin lines represent the resulting functions when the 
photometric error was assumed to be e=0-0^ mag in images A/B and (ts£;=0.02 in images C/D, and (t^s=0.02 mag in images A/B 
and (ri^B=0.04 in images C/D. The lens model contained 50% smoothly distributed matter. 




Figure 4. The functions representing the probability for combined image event rates in Q2237+0305. Left: the time averaged rate, with 
expectation values, Centre: events expected per decade and Right: events expected per 6 month observing season. The light and dark 
lines refer to assumed directions for transverse motion that are directed along the A-B axis and C-D axes respectively. The thick and 
thin lines represent the resulting functions when the photometric error was assumed to be (tsb=0.01 mag in images A/B and (Tg£;=0.02 
in images C/D, and ai^E=0.02 mag in images A/B and crLE=0.04 in images C/D. The lens model contained 50% smoothly distributed 
matter. 



bution. There is circumstantial evidence for approximately 
5 HMEs in existing light-curves for Q2237-I-0305 (Irwin et 
al. 1989; Corrigan et al. 1991; 0stensen et al. 1995; OGLE 
Web page). This number is larger than, but consistent with 
the mode of most of the distributions. However, it may be 
improbably high if the combination of a trajectory along the 
A-B axis and 50% smooth matter content is assumed (for 
both assumptions of photometric error). 

There is 1 chance in ~ 4 — 10 (20 for cases where a tra- 
jectory along the A-B axis and 50% smooth matter content 
are assumed) that a SHME will be observed during a given 
observing season. However the observation of two SHMEs 



in one observing season is very unlikely (< 1/25), indicat- 
ing that the temporal scale of a diamond shaped caustic 
is typically larger than 6 months. Witt, Kayser & Refsdal 
(1993) found that M-shaped double events have an aver- 
age separation of 0.9-1.3 years (they assumed a source plane 
transverse velocity of 6000 fcm see" ^ and a mean microlens 
mass of ~ 0.35). This is demonstrated by the probable dou- 
ble peaked event in image A (Irwin et al. 1989; Corrigan et 
al. 1991; Racine 1992) which has a duration of ~ 500 days. 

Figure ^ shows histograms of the number of SHMEs ex- 
pected in the 10 years immediately following the observation 
of an HME. The figure shows a large excess in the rate of 
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Figure 5. The functions representing the probability for events expected in the decade following the observation of an HME. The light 
and dark lines refer to assumed directions for transverse motion that are directed along the A-B axis and C-D axes respectively. The 
thick and thin lines represent the resulting functions when the photometric error was assumed to be CTg£;=0.01 mag in images A/B and 
asE=0.02 in images C/D, and o'/,£;=0.02 mag in images A/B and (t^£;=0.04 in images C/D. 



events during such periods. Indeed, the most likely number 
of events is non-zero. The means of these histograms, as well 
as the correspondine means computed for the average rate 
are shown in table y. The values demonstrate a ~ 50 — 100% 
increase in SHME rate in the decade following an HME over 
the average value. 

Two images in Q2237+0305 have undergone a signifi- 
cant change in brightness between the end of the 0stensen et 
al. (1996) data and the beginning of the 1998 observing sea- 
son. The expected rate of SHMEs for the next 10 years may 
therefore be significantly higher than the values presented 
in table eI 



5 CONCLUSION 

We have combined calculations of the time averaged rates 
of SHMEs over a phase space of assumed values for effec- 
tive galactic transverse velocity, mean microlens mass and 
continuum source size with previously obtained probabili- 
ties for their values. We find that the total SHME rate in 
Q2237+0305 is 1.5±0.6 - 6.2±1.3 events per decade depend- 
ing on the assumptions made for the direction of the galactic 
transverse velocity, the smooth matter content and photo- 
metric error. We find that the probable number of observed 
HMEs is consistent with calculated distributions of events 
per 10 year period. During periods following a caustic cross- 
ing in an image, we find that the event rate in that image 
is between ~50 and 100% higher than the average value. In 
addition, we have calculated the distribution of SHMEs that 
should be observed per 6 month period and find that there 
is one chance in ^ 4 — 10 that a SHME will be observed. The 

6 month period corresponds approximately to one observ- 
ing season. These results therefore add further support to 
the case for continued monitoring of Q2237-I-0305 with the 
aim of making spectroscopic and photometric observations 
of Q2237-I-0305 during a SHME. 
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